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A rigorous solution for the spectrum of a quasioptical cylindrical cavity resonator with a randomly 
rough side boundary has been obtained. To accomplish this task, we have developed a method for 
the separation of variables in a wave equation, which enables one, in principle, to rigorously examine 
any limiting case — from negligibly weak to arbitrarily strong disorder at the resonator boundary. 
It is shown that the effect of disorder-induced scattering can be properly described in terms of two 
geometric potentials, specifically, the "amplitude" and the "gradient" potentials, which appear in 
wave equations in the course of conformal smoothing of the resonator boundaries. The scattering 
resulting from the gradient potential appears to be dominant, and its impact on the whole spectrum 
is governed by the unique sharpness parameter S, the mean tangent of the asperity slope. As 
opposed to the resonator with bulk disorder, the distribution of nearest-neighbor spacings (NNS) in 
the rough-resonator spectrum acquires Wigner-like features only when the governing wave operator 
loses its unitarity, i. e., with the availability in the system of either openness or dissipation channels. 
It is shown that the reason for this is that the spectral line broadening related to the oscillatory 
mode scattering due to random inhomogeneities is proportional to the dissipation rate. Our numeric 
experiments suggest that in the absence of dissipation loss the randomly rough resonator spectrum 
is always regular, whatever the degree of roughness. Yet, the spectrum structure is quite different 
in the domains of small and large values of the parameter 3. For the dissipation-free resonator, 
the NNS distribution changes its form with growing the asperity sharpness from Poissonian-like 
distribution in the limit of H ^ 1 to the bell-shaped distribution in the domain where S > 1. 

PACS numbers: 42.60.Da, 32.30.Bv, 68.49.-h 



I. INTRODUCTION 

A great deal of attention is currently being given to 
studying the spectral properties of microwave systems in 
which the propagation of waves is highly, or even totally, 
restricted in all three dimensions. Systems of this sort, 
widely referred to as zero-dimensional systems, are com- 
monly used in optic, acoustic, and laser technologies as 
resonators, both of cavity and of open type. 

Spectral properties of resonators whose spatial dimen- 
sion substantially exceed an operational wavelength (the 
so-called quasioptical resonators) have for a long time 
been thoroughly investigated (see, e.g., Refs. ||). 
However, this applies mainly to resonators of relatively 
simple form, which allow for exact intcgrability of corre- 
sponding wave problems. Things are getting worse with 
resonators whose geometrical features do not permit the 
solution of wave equations using conventional variable 
separation methods. Such resonators, provided they have 
no internal inhomogeneities, arc, by analogy with classi- 
cal dynamic systems, attributed to the class of nonintc- 
grable systems of the billiard type. 

In studies of billiard dynamical systems, the most com- 
monly used are the methods of classical chaos theory [|| , 



'Electronic address: yutarasov@ire.kharkov.ua 



which rely substantially on the theory of random matri- 
ces |Q, Using these methods, spectral properties of 
different billiard systems with unstable classical dynam- 
ics have been profoundly studied, in particular, the well- 
known Sinai and Bunimovich billiards [?j . Systems of 
similar type, with regular and mainly smooth boundaries, 
were examined experimentally using different microwave 
methods ||. At the same time, spectral properties of 
resonators whose boundaries possess essentially random 
properties have thus far received insufficient attention. 
The adequate theoretic analysis of these properties is 
difficult to perform, mainly, due to the widespread be- 
lief that in systems with broken spatial symmetries the 
variables do not separate, at least with a traditional alge- 
braic approach. This makes it necessary to employ rather 
approximate, and even phenomenological, methods for 
spectral analysis of such systems, whose application is 
far from being adequately substantiated. 

In our recent papers, Refs. |^-p"l[, using the method 
developed previously in Ref . |12 f we have shown that one 
can separate variables in a wave equation and thus per- 
form a comprehensive spectral analysis, for any restricted 
wave system, however disordered, lossy, or lossless it may 
be. The only requirement for the system in question is 
that its difference from some deliberately integrable ref- 
erence system might be described in terms of certain ef- 
fective potentials. In this case, the variables can always 
be separated, although not necessarily by the algebraic 
approach but rather with some operator technique. Nev- 
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crthclcss, this enables one to scrutinize any limiting case, 
from negligibly weak to arbitrary strong scattering of os- 
cillatory modes due to inhomogencities in the resonator. 

In the above-mentioned papers it was found that the 
scattering produced by static inhomogeneities smoothly 
distributed in the resonator bulk, though Raylcigh in na- 
ture, results in additional spectral line broadening, which 
is of essentially nondissipative origin and is additive to 
the broadening caused by dissipation processes. This 
kind of broadening is explained by the resonance mode 
dccohcrence resulting from coexistence in such systems of 
two scattering channels, specifically, the intramode (co- 
herent) and the intermode (incoherent) channel. The im- 
pact of infill inhomogeneities on the separate resonance 
lines is of a highly selective character. The most widened 
and shifted along the frequency axis appear to be those 
lines in whose vicinity there is the largest number of other 
resonances. Meanwhile, the lines positioned relatively 
sparsely prove to be weakly subjected to inhomogeneities 
and retain their high quality factors. 

Such a specific behavior of spectral lines was inter- 
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as a kind of interaction between 



preted in Rcfs. 

resonances. The interaction results in effective rarefica- 
tion of the initially dense spectrum of the random inho- 
mogeneous resonator, in the sense that as the disorder 
increases the number of high-Q-factor lines gets progres- 
sively lower. The quality factors of the remaining lines 
are considerably reduced, and some of the lines even be- 
come superposed. We have exploited such an effective rar- 
cfication of cavity resonator spectrum, which results from 
randomization of its infill, in simulations of a random 
laser with single-frequency generation using the quasiop- 
tical millimeter wave band resonator Jf3[ . 

In this paper, the investigation is undertaken of the 
cylindrical quasioptical cavity resonator which is ran- 
domized not through inhomogeneous filling but rather 
through random roughening of its side boundary. Res- 
onators of such a type are physically equivalent to 
the two-dimensional billiard system conventionally called 
the Shepclyansky billiard. It was already examined in 
Refs. I4-M in an effort to elucidate the possibility to 
observe the effects of dynamic localization and quantum 
ergodicity. Apart from this, Refs. |l^, [l9| also dealt with 
resonance systems having randomly rough boundaries, 
specifically, with solid-state quantum dots. To study the 
current-carrier spectrum, the authors of Refs. jl8|, [l9tl 
have applied the ballistic a model intended to describe 
not only universal but also system-specific properties of 
quantum resonance systems. Yet, carrier scattering due 
to boundary roughness was taken into account in these 
papers phenomenologically, using the concept of the spec- 
ular reflection factor |2(|. This made it impossible for 
the authors to allow for the diffraction effects produced 
in the course of quantum wave scattering and, as a con- 
sequence, to carry out in-depth spectral analysis for the 
system where the boundary scattering is of determinative 
significance. 

Some more attempts of theoretical description of 



rough- wall resonator spectra can be noted. However, 
they all run inevitably into great difficulties associated 
with the lack of appropriate methods for considering 
the wave scattering at random inhomogeneous bound- 
aries of the system under investigation. For instance, 
in spectral analysis of open-type resonators possessing 
randomly rough boundaries, namely, dielectric disk res- 
onators (DDRs), the method of polarization currents (the 
so-called volume current method, VCM) was widely uti- 
lized in due time [ pl| . Primarily, the method was in- 
vented to make plausible estimations of radiation loss 
of whispering gallery (WG) oscillations of randomly 
bounded open electrodynamic systems [^2], ^3|. This 
method reduces essentially to simulation of oscillation 
scattering produced by real boundary inhomogeneities 
by adding the fictitious sources (currents) randomly dis- 
tributed near the unperturbed (nonrough) boundary of 
the system. However, in our recent paper, Ref. pi) , 
we put forward the arguments in favor of the inconsis- 
tency of using the VCM to describe the radiation loss in 
DDRs with rough boundaries and obtained the solution 
of this intricate problem with a more efficient and reliable 
method, namely, the direct solution of wave equations 
for fields in such systems. This was made possible due 
to elaboration of the specific method for resonance mode 
separation, which is applicable to systems with arbitrary 
degrees of disorder. We succeeded in showing that the de- 
cay in quality factors of resonators with randomly rough 
boundaries results, basically, from the roughness-induced 
intermode scattering. Due to this type of scattering, the 
energy of high-Q WG modes goes over to other, less sta- 
ble, modes and in such a way radiates out of the resonator 
volume. Interestingly, the level of radiation loss of a ran- 
domly rough DDR turns out to be mainly governed by 
the average slope of the asperities (or, in other words, by 
their sharpness) and is much less sensitive to the ampli- 
tude of boundary fluctuations. 

It should be noted that taking into account the random 
roughness of wave system boundaries is a rather complex 
mathematical problem, which has by now a long-lasting 
history. Basically, different methods for describing the 
wave scattering by such boundaries were previously de- 
veloped for open and infinite systems, such as sea sur- 
face p5|-p7[. Even so, the number of problems subject 
to being completely solved was principally restricted by 
the applicability of the small-perturbation method (small 
amplitude and small slope of boundary asperities) and/or 
Kirchhoff approximation (the inclusion of one-fold scat- 
tering in the geometrical-optic approach). As to the essen- 
tially confined systems such as resonators and quantum 
dots, their spatial boundedness imposes, on its own ac- 
count, significant limitations on the possibility to apply 
statistical methods for describing the wave scattering in 
such systems. Therefore, having in mind the essential 
finitcness of the object of our study in this paper, the 
methods used previously to cover wave scattering caused 
by boundary roughness need serious modifications. 

To calculate the spectra of resonance systems with ran- 
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domly inhomogeneous boundaries we suggest a method 
which makes it possible, above all, to perform calcula- 
tions without imposing the substantial limitations on the 
random asperity amplitude. Second, the method enables 
one to avoid any restrictions regarding the asperity sharp- 
ness. The particular advantage of the suggested method 
is that it is not rigidly bounded to symmetries of the 
system under consideration. 

The method includes two basic stages. Initially, the 
problem regarding the oscillations in cylindrical res- 
onator with arbitrarily rough walls is reduced, by means 
of the boundary conformal smoothing, to the problem for 
oscillations in the resonator whose boundaries are ideally 
homogeneous. 1 48 In this way, one is able to introduce a 
new complete set of eigenmodcs (we call them gradient 
renormalizcd, or GR, modes) which serve as a good basic 
approximation for developing the appropriate perturba- 
tion theory. 

At the next stage, the GR modes are entirely sepa- 
rated, much in the same way as it is done for standard 
integrable systems. Remarkably, however, these specific 
modes can be separated regardless of the level of the as- 
perity sharpness, the latter being the measure of scat- 
tering degree. The deduced functional expressions for 
cigcnfrequencies of the disordered resonator enable one 
to carry out a detailed analysis of its spectral properties 
to any desirable accuracy. 

The conclusions we have made from our considerations 
are as follows. First of all, at the amplitude of bound- 
ary inhomogeneities, which is small as compared to res- 
onator dimensions, the level of the disorder is specified 
by the universal gradient parameter, namely, the average 
slope of the asperities, rather than by the relationships 
between the mean asperity amplitude, the asperity cor- 
relation length, and the wavelength of the excited oscil- 
lations. Depending on whether the asperities are smooth 
(i.e., their mean slope ratio is small as compared to unity) 
or sharp, the resonator behaves either as a weakly or a 
strongly disordered wave system. 

Next, as opposed to resonators with the disordered 
bulk infill, the distribution of nearest-neighbor spacings 
in the spectra of surface-inhomogeneous resonators ac- 
quires the features typical for Wigncr distribution (which 
is conventionally associated with chaotic properties) only 
when the unitarity of wave operator is violated, i.e., if 
dissipation and/or openness channels arc available in the 
system. In the absence of dissipation the spectrum is 
completely regular, with zero-width resonance lines, but 
offers essentially different correlation properties and reso- 
nance level densities in the domains of smooth and sharp 
boundary asperities. 

The result of particular interest we have achieved in 
the present work is the relationship established between 
the roughness-induced decoherence of resonance modes 
and the resonator dissipative properties. The resonance 
line broadening related to (Raleigh-type) mode scatter- 
ing due to randomly rough boundaries, being additive 
with the broadening caused by the dissipation per se, is 



found to be also proportional to the dissipation rate. If 
dissipation (or resonator openness) is nonexistent, so will 
the roughness-induced broadening be unavailable. This 
fact fundamentally differentiates resonators with bound- 
ary inhomogeneities, which belong to the class of billiard- 
type systems, from resonators randomly inhomogeneous 
in bulk. 



II. THE CHOICE OF THE MODEL AND THE 
PROBLEM FORMULATION 

The object of our study is a resonator of cylindrical 
form with generatrix parallel to z-axis (see Fig. p, which 
is bounded by ideally-conducting plates passing perpen- 
dicular to z-axis through the points z± = ±H/2 and the 
ideally-conducting side boundary 



S = 



r = R + Z(V>) 

if G [ — 7T,7r] 



(1) 



Here R is the average radius of the resonator, i.e., the 
radius of an ideally circular cylinder whose volume co- 
incides with that of the original rough cylinder having 
the same height. The function £(</?), which specifies in- 
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FIG. 1: Schematic view of cylindrical resonator with a ran- 
domly rough side boundary. Vectors n and no point toward 
outer normal to the rough and ideal boundary, respectively. 

homogeneity of the resonator side boundary, is assumed 
to be the random process with zero mean, (£,(<p)) = 0, 
and binary correlation function 



{(L{<pW)) = ° 2 w{v-v') 



(2) 



Taking a in Eq. (g) to be the mean-square height of the 
roughness, correlation function W(ip) will be thought of 
as having unity maximum at <p = and decreasing to 
paramctrically small values on the angle scale | A^l = ip c . 
The correlation angle <p c is related to the roughness arc 
correlation length, s c , through the obvious relationship 
(p c = s c /R. Both of the functions, (,(ip) and W(ip), are 
taken to be periodic with the period 2n. In what follows, 
we regard the roughness as small scale in that the arc 
correlation dimension is small as compared to the char- 
acteristic size of the system, 



< R 



(3) 
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Despite the resonator being a restricted system, this 
inequality suggests the conditions for spectrum self- 
averaging to be fulfilled with parametric accuracy. An 
additional argument in favor of the self-averaging is the 
well-known fact that classical dynamics in systems anal- 
ogous to wave billiard under study, whose boundary per- 
tains to the class of scattering boundaries, is essentially 
ergodic ||, |0 . 

Upon finding the Greens function G(r, r') of the 
Helmholtz equation inside the region depicted in Fig. [I], 
by way of boundary conditions (BC) we use the equalities 



1 



3G(r,r') 



dr„ 



and 



G(r,r') 



= , 



(4a) 



(4b) 



z=±H/2 



which correspond to .ff-polarized electromagnetic oscil- 
lations (TE oscillations). Q The derivative in Eq. ( |ia| ) 
is taken over an outer normal line to actual, corrugated 
side boundary (see Fig. 111). 



In the general case the solution to wave equations with 
BC Eqs. (||) is a challenging task. In this paper, following 
the approach we have previously applied in Ref. ]24| , we 
smooth out the rough boundary of the resonator to the 
ideal one, i.e., unperturbed by the asperities. In this 
way we manage to reduce the problem of wave scattering 
at the resonator actual boundary to the scattering by 
fluctuations of the effective medium filling the bulk of a 
fictitious resonator, the cross section of which is the ideal 
circle. 

We relate the coordinates in the initial cylinder system 
coupled to the resonator axis to new cylindric coordinates 
(we mark them by the tilde sign) through the conformal 
transformation 



Y 

Z 



v > 



(5) 



whose Jacobian is equal to 1/ [1 + £((p)/R]. In the new 
coordinate system the Green function equation assumes 
the form 



ld_ d_ 

r dr dr 



i d 2 d 2 



dip 2 dz 



J 



+ — + K 2 - V {h) - V {s) G{r, (p, z; r'tp', z') = -5(r - r')5{ip - <p')5(z - z') 



(G) 



(the tilde signs are omitted from here on) , differing from 
the initial equation by the additional terms in the wave 
operator, namely, V^ h > and V^ s \ which we call below the 
effective potentials, 



k 



dz 2 
d 



i?/%) dip 



g §M ' 



ld_ 

r dr 

n 2 



ld_ d_ 

r dr dr 



(7a) 



(7b) 



The analytical structure of these potentials is determined 
by the roughness function itself, £(ip), and by its deriva- 
tive £'(<£>). In Eq. (||), the notation K 2 = k 2 — i/T dis 
is introduced, where k — ut/c and 1/Vdi s is the effec- 
tive dissipation rate which we put in phenomcnologically 
[time dependence of the fields is chosen to be of the form 
exp(— iuxt)]. Factor ft(ip) in Eqs. (|^) is used to simplify 
the formula, 



1 



R 



(8) 



Upper indices "ft" and "s" at the potentials (|7|) symbolize 
the fact that one of them, namely, the former, is governed 



by fluctuations of the asperity height solely [i.e., directly 
by function £(<£>)] whereas the potential is mainly 
determined by the asperity slope [by the derivative £'('y 9 )]- 
We consider £(y>) to have a relatively small statistical 
value and put /3(<p) ~ 1 below. 

In view of inequality (||), one can neglect the correla- 
tion between functions £(y>) and at angle distances 
larger than tp c . Owing to this we can relate two mu- 
tually independent physical mechanisms with potentials, 
Eqs. (R), which are responsible for wave scattering by sur- 
face roughness, namely, the "height" and the "gradient" 
mechanism, both predicted previously in Ref. plj |. 

Note that contrary to the common belief (see, e.g., 
Ref. |2j| ) , even though we hold the asperities to be small 
in height in the sense of inequality 



cr < R 



(9) 



the effective potential (^) may take on quite large val- 
ues. The local value of this potential and, consequently, 
of the entire potential in Eq. (^), is determined by ra- 
tio £'(fp)/R, which under constraint (0) assumes arbi- 
trary large values. Clearly, this does not necessarily sig- 
nify large scattering strength, as the latter is the integral 
rather than the local property of the potential. Below we 
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formulate the particular parametric conditions enabling 
one to regard the effect of scattering due to effective po- 
tentials (pi) as either weak or strong. 



III. GRADIENT RENORMALIZATION OF 
WAVE OPERATOR AND SEPARATION OF 
OSCILLATION MODES 

In contrast to "hight" potential V^ h \ "slope" potential 
V( s ) in Eq. ([}]), keeping in mind the prospective usage of 
perturbation theory, is defined in an awkward fashion 
as its average is not equal to zero. By separating this 



average as the operator summand, 

r or or 



where we have introduced the gradient parameter 



(10) 



(11) 



which specifies the asperity sharpness degree, we can re- 
cast equation (ra) in the form 



J 



„ „ 2 ,1 d d 1 d 2 
(1 + ^ 2N 



if 



r dr dr r 2 dtp 2 dz 2 



+ K 2 - V {h) - V {sl) - V" (s2) 

r 



G(r,<p,z;r'<p',z') = -S(r-r')S(ip-<p')5(z-z') . (12) 



In Eq. (|l2|), two different slope potentials are introduced 
instead of the potential (7b), whose configurational aver- 
ages arc by definition equal to zero, namely, 



1 



R 

y(s2) = _ 



ld_ 

r dr ' 



ld_ d_ 

r dr dr 



(13a) 
(13b) 



The function P(<p), subject to proviso (Q), is set equal to 
unity herein. 

Let us proceed to Fourier representation using thereto 
a complete set of eigenfunctions of the differential op- 
erator in Eq. (|l^), which are consistent with boundary 
conditions (Q) and (^-periodicity of the desired solution. 
In Eq. (|l^) with no random potentials, the variables are 
readily separated, thereby giving the following complete 
set of normalized eigenfunctions: 



|r ;/ x> = |^n>|r;0 (S) M 
(n = n/y/l + S 2 ) ; 

\<p;n) = (2 7 r)- 1 / 2 e mv 
(n = 0,±l,±2, ...) , 

| r;Z) (n) =jB (n)V| J| _ |(7 p|) r/jR) 



(1 = 1,2, ...) 



\z;q) = W^sin 
(q=l,2, ...) 



Here, 7 ; ""'' 1 denote the positive zeros of function </|~|(i) 
enumerated by natural numbers (I) in ascending or- 
der. The normalization constant in "radial" eigenfunc- 



(14a) 

(14b) 
(14c) 

(14d) 



tion (14c) is 



D 



an 



Ji 



(l»l) 



Jl 



(N) 



1 -1/2 



(15) 

Below we refer to the Hcrmitian differential operator 
m Eq. (P), which possesses the set of eigenfunctions 
(p~4]), as the GR Laplace operator. In representation of 
this operator eigenmodes Eq. ( fl2"|) reduces to a set of 
non-homogeneous coupled equations against the Green's 
function Fourier components, 

(K 2 V M ) - U^G V ^ = 5^ . (16) 



Here. 



G. 



drdr'(r;n\G(r,v')\r';n') (17) 



is the Green's function in momentum representation, and 
il is the domain in coordinate space M 3 occupied by the 
resonator with an unperturbed boundary. The quantity 
>f 2 , which is the GR Laplasian eigenvalue, in what fol- 
lows is called an unperturbed energy of the /xth resonator 
mode, 



^ = (l + 5 2 



7* 



(N)' 



R 



(18) 



Accordingly, the eigenfunctions ([14]) will be symbolically 
termed as unperturbed eigenfunctions, even though they 
contain gradient parameter 5. 

Matrix of the Eq. ( |l6| ) coefficients is composed 

of matrix elements of the total perturbation potential in 
Eq. (|l2"|), which contains both the height potential and 
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the gradient ones [V(P ert ~>(r) = + + V^ s % 



U^= f dr(r; / i|V' {pert) (r)|r;i/> 



(19) 



The diagonal element of this matrix, = V M , is de- 
liberately separated from the sum in Eq. @. The scat- 
tering induced by this element is specified below as the 
intramode, or coherent, scattering, and the potential 
will be named the intramode potential. On the contrary, 
all the constituents of matrix IIZ^^H with distinct mode 
indices (/i ^ v) are thought of as responsible for incoher- 
ent scattering between different modes and are nominally 
referred to as intermode potentials. In Refs. ]l2], |24|, [3^] 
it was shown that such a segregation of the whole set 
of potentials into "coherent" and "incoherent" subsets 
enables one to correctly deduce the whole set of uncou- 
pled equations for diagonal elements of the entire Green's 
matrix ||G,j,„||, through which all other elements are then 
expressed by linear operation. Basically, this fact makes 
it easy to analyze the quantum (wave) systems subjected 
to quite arbitrary external potentials, including nonlocal 
ones. 

Equation ( |l6|) allows the peculiar separation of vari- 
ables provided that at least one of two fundamental sym- 
metries, either V or T, is broken. The variables are 
separated through the operator procedure detailed in 
Ref. jllj. The key point of the procedure is to intro- 
duce a certain trial mode propagator G^ for each of 
the harmonics, say, with mode index u. This propagator 
represents the "seed" mode Green's function, i. e. the so- 
lution to Eq. ( |l6| ) obtained without regard to intermode 
potentials, and, hence, intermode scattering, 



(20) 



Starting from this trial function, all the intermode prop- 
agators G vfJl (with \/u ^ n) are expressed in terms of 
just one intramode Green's function by means of 

functional equality: 



G V fj, = P„{t - R) 



1 R-F^G^ 



(21) 



Here, operator R has the meaning of a mode-mixing op- 
erator, whose domain of definition is mode subspacc 
containing the whole set of mode indices save index fi; P^ 
is the projection operator that assigns the given value fj, 
to the nearest mode index of any operator standing next 
to it, no matter if it is to the left or to the right. The 
operator R is expressed in the product form, 



r = g^u 



(22) 



where operators and U are specified on M M by their 
matrix elements 



{v\U\u') = U V u> ■ 



(23a) 
(23b) 



Putting then mode index fi' = \i in Eq. ( |l6| ) and sub- 
stituting all the intermode propagators in the form (21), 
we arrive at the infinite set of uncoupled equations for all 
intramode propagators G^, whose solution is 



Guv 
Here, 



[K 2 



(for V/Vs) . (24) 



T M = P fJ ,U(l-R)- 1 RP fl 



(25) 



is the portion of the /xth mode eigenenergy which is re- 
lated to the intermode scattering only. The whole set 
of intramode propagators, Eq. (^4|), in conjunction with 
relationship Eq. (Ell), solves entirely the perturbed res- 
onator Greens function. In such a way we have re- 
duced the problem of determining the randomly rough 
resonator spectrum to finding the poles of solely diago- 
nal elements of matrix |jG Miy ||. 

From Eq. (^j) it can be easily seen that in the gen- 
eral case all eigenfrequencies of the resonator we study 
appear to be complex valued. The complexity is appar- 
ently introduced by two sources. One of them is related 
to dissipation directly, being accounted for through the 
dissipation rate entering complex energy of the resonator, 
K 2 = k 2 — i/Tdi S . The other originates from the T po- 
tential (|2^), which is nothing but the diagonal (in mode 
representation) element of the regularized T matrix well- 
known in quantum scattering theory p3| , p4[ . This po- 
tential functional structure suggests that it can be in- 
terpreted as taking into account the effective interaction 
between the particular mode fi and the rest of trial modes 
with Vf ^ fi. 

Eventually, with the account of intermode interaction, 
the mode fi becomes either coherent or incoherent, de- 
pending on whether the dissipation is present in the res- 
onator (1/rdis 7^ 0) and/or the potential ( |25| ) contains 
an imaginary part. The latter kind of imaginarity comes 
about in the self-energy of the particular wave or quan- 
tum system (in our case, the mode state /x) interacting 
with the hardly controllable "environment" if one tries 
to describe this system separately. Normally the imagi- 
nary addendum to the isolated system self energy is in- 
terpreted as environment-induced indefiniteness of its en- 
ergy levels, or, in other words, as the result of its deco- 
herencc j3f| (see also Ref. [^6|). Hereafter we also refer 
to the imaginary part (if any) of the T potential ( p5| ) as 
the decoherence rate of the /xth eigenlevel. 



IV. STATISTICAL ANALYSIS OF THE 
RANDOMLY ROUGH RESONATOR SPECTRUM 

Since the steady-state propagation of waves in the res- 
onator under study is expected to be ergodic (see Ref. 1 37 
and references therein), one can obtain its spectrum by 
performing configurational averaging of the Green's func- 
tion (pif). We denote this type of averaging by angular 
brackets, (■)• 
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The net random potential in Eq. (Q) and, consequently, 
its matrix elements specifying the intra- and intermodc 
potentials in Eq. ( pi| ) , are strongly dependent on the sta- 
tistical properties of the roughness function £((f). The 
intramode scattering produced by the potential V M may 
be regarded as either weak or strong, subject to small or 
not, in comparison with the unperturbed mode energy 
>C, appears to be the mean-square norm of this poten- 
tial estimated on the unperturbed mode eigenfunction. 

As far as intcrmode scattering is concerned, the asso- 
ciated additions to the unperturbed spectrum, Eq. (fl8|), 
result mainly from the potential 7^, whose average is not 
equal to zero. This potential has a rather involved func- 
tional structure; therefore it seems to be reasonable to 
characterize its effect by estimating the norm of operator 
R which controls the trial mode intermixing. Note that 
these modes are not the eigenstates of a parent nonrough 
resonator. Their energies, as seen from Eq. (^|), contain 
the sharpness parameter S 2 , and, therefore, the field of 
each of these modes may be represented as a bundle of 
parent, nonrough resonator modes. Keeping this fact in 
mind, we call the above-introduced trial modes, which are 
specified by eigenenergies (|l8|), the unperturbed compos- 
ite modes. 



A. The comparative estimation of scattering 
mechanisms for composite modes 

1. Intramode scattering 

Under restriction (^|), the height potential, Eq. (|7a|), 
can be represented with good accuracy by approximate 
expression 



R 



(26) 



Since the random process £(<p) is supposed to be ergodic, 
the diagonal ma trix element of potential (|2(]), taken be- 
tween functions ( |l4rj| ) and its complex conjugate, is equal 
to zero. Diagonal matrix elements of either of gradient 
potentials ( |l3| ) also vanish: (i) for the potential V^ s1 ^ 
this is because of its periodicity in angle ip, and (ii) as 
regards the potential V^ s2 ^ this is due to the assumed er- 
godicity of the roughness function. |30[ Hence, under our 
assumptions regarding the statistical properties of 
there is no coherent scattering of composite modes which 
would be caused by random roughness of the resonator 
side boundary. 



2. Intermode scattering 



functional structure of this matrix and, consequently, the 
magnitude and the analytical structure of the potential 
Tfj, are determined by the value of the mode-mixing op- 
erator R. We regard the intermode scattering as either 
weak or strong depending on whether the average norm 
of this operator is small or large as compared with unity. 
In view of the operator R multiplicative structure, we will 
estimate this norm using the formula 



<||R|| 2 ) = max ( I £GMW„ 



(27) 



Here, trial Green function Gv is independent of the 
roughness function with regard to remarks made in 
the previous subsection. The averaging of the potential 
U vlil modulus square, which is in fact simple but requires 
a gret deal of tedious calculations, is detailed in the Ap- 
pendix. Eventually, we obtain the following estimations: 



<|jRW|| 2 )~(a/i?) 2 , 
<||R^!| 2 )~(a/ Sc ) 2 ~~ 2 
<||R( 82 >|| 2 >~(<7/i?) 4 . 



(28a) 
(28b) 
(28c) 



From these estimates one can infer that, in view of 



inequality (|9|), the potentials U$ and U^v' result in 
negligibly weak intermode scattering. For the scattering 

induced by the potential U$ , it can be easily seen 
that under conditions (|^) and (^|), first, the norm of 
the corresponding term in the intermixing operator, 
R( sl ), is the uppermost among the norms presented in 
Eqs. (p8|), and second, the operator R( sl ) norm can take 
on both small and large values, as compared to unity. 
This enables us to analytically treat any of the limiting 
cases, both of weak and of strong intermode scattering. 



As2) 



B. Evaluation of spectrum corrections 



The composite mode spectrum given by Eq. (|18| ) , even 
without additional corrections induced by potentials en- 
tering Eq. ( |l6|) , can differ essentially from the initial, un- 
perturbed resonator spectrum. The difference is either 
small or large depending on whether the mean asper- 
ity slope (in fact, the parameter S 2 ) is small or larg e 
as compared to unity. As may be seen from Eq. (28b), 
it is exactly the same parameter that governs the level 
of composite mode mixing. Below, while analyzing the 
spectral corrections, we keep in mind only the ones that 
result from the scattering induced by the potential V^ sl \ 



The effect of intermode scattering on the resonator 
spectrum, that is, on the pole structure of function (pi[), 
is governed by the potential (|5j), which is defined as the 
diagonal matrix element of the regularized T matrix. The 



1. Weak scattering limit 

As it follows from estimates (p8|), the intermode scat- 
tering, or, in other words, the entanglement of composite 
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resonant modes, is weak provided that two inequalities 
hold simultaneously, namely, 



a/R < 1 
<t/s c < 1 



(29) 



meaning the smallness (on average) of both the ampli- 
tude and the slope of random asperities. In this limiting 
case, one can expand the inverse operator in Eq. ( |25| ) 
in power series in the operator R, thereby obtaining the 
approximate expression 

T M « PjAG^UP^ = U^G^U^ (30) 

for the diagonal element of the T-matrix. The correlator 
of intermode potentials in Eq. (|o|), subject to definition 
(||) and estimates (p8|), is calculated as 

x W( ntl - n v ) (nl - nlfv^V,,, , (31) 



where 2? M „ is the matrix element of the opera- 
tor ( 1/r) (d/dr), which is taken between eigenfunc- 
tions (14c), 



TV 



dr{r;l„\^-^\r;l v )^ . (32) 



By representing the average mass operator for mode 
fi as a sum of real and imaginary parts, (7^) = 

Ax^ + i'/r$ , we arrive at the following expressions for 
resonance line shifting and broadening, which are associ- 
ated with composite mode scattering, 



2 2 

< - K 



^-^) 2 + (iAdi S ) 2 ' 

l/Tdis 



The weight factor in Eqs. ( |33| ) is given by 



(33a) 
(33b) 




Interestingly, this factor is not positive definite in the 
general case. If we consider, instead of TE polarized oscil- 
lations, the oscillations of TM type, for which boundary 
conditions arc applicable 



G(r,r')| 







and 



9G(r,r') 



(35) 



z=±H/2 



rather than BC Eq. (g) (see Ref. f§), then the weight 
factors appropriate to the transitions between any two 
different modes would be necessarily positive. This stems 
from the fact that for TM oscillations the quantities 7}'"" 
in Eq. (|34|) would be zeros of the Bessel function itself 
rather than zeros of its derivative. Based on this fact, 
antisymmetry of the tensor T)^ would follow and, corre- 
spondingly, the antisymmetry of each of the integrals in 
Eq. (p4j). As for TE oscillations being considered here, 



the sign of the coefficient for different pairs of mode 
indices is not a priori definite. This implies that the in- 
termode scattering contribution to Q factors and shifting 
of different lines can be of different sign. To verify this 
anticipation by analytical means seems to be an unreal- 
istic task. Therefore we subsequently analyze the rough- 
resonator spectrum numerically, using for this purpose 
Eqs. (||) and @. 

Yet, even at this stage we are in a position to make 
some qualitative statements regarding the influence of 
surface inhomogeneities on the resonator Q factor as well 
as on the position of resonance lines. Assuming the az- 
imuth mode indices to be large as compared to unity, 
one can replace exact zeros 7; with their asymptotic 
expre ssion s (see Ref. @). For scattering-induced decay 
rate (|33b|), the approximative formula readily follows 



—(- 



[R 2 /T di „ 



(36) 
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Here, the sum over radial mode indices, subject to the re- 
striction imposed on the summation region in Eqs. (j33|), 
is evaluated to be of order 1 + (i? 2 /rdi s ). The conver- 
gence of the sum over azimuth indices is provided by the 
correlation function ^(n^ — n„), from which we obtain 



1 

Tdii 



S 2 (n, 



_-n 2 



(37) 



Because each of the resonance line Q factors is deter- 
mined by the sum of the dissipation rate and the rate of 
roughness-induced attenuation, 



o; 1 



^l/r dis + 1/t< sc) 



(38) 



the emergence of random inhomogeneities on the res- 
onator side boundary may result either in the preser- 
vation of the specific line width, if the condition 
S 2 (n^ + tfic 1 ) 1 holds true, or in its significant broad- 
ening, with this inequality changing to the opposite one. 
It is, however, intriguing that Q factors of all resonance 
lines, subject to inequality (^|), remain infinite in the ab- 
sence of dissipation, irrespective of whether the resonator 
side boundary is rough or not. At the same time, the po- 
sitioning of these lines can differ considerably from that 
pertinent to the initial, nonrough resonator. Even if the 
asperities may be very smooth (S <C 1), the shift of the 
specific line, being determined by the difference between 
energy (|l8| ) and the analogous energy calculated for 3 
exactly equal to zero, can substantially exceed the mean 
level spacing. In principle, this single fact suffices to 
change significantly the distribution of resonances along 
the frequency axis, as compared to the inhomogeneity- 
free resonator. But apart from that, the additional shift 
should also be taken into account, which is governed by 
Eq. ( |33aD , whose relative value, as numerical estimations 
do reveal, may significantly exceed the parameter S 2 in 
view of a large number of terms in the sum. 

Of dominant importance, however, is the fact that the 
spectrum of weakly rough [in the sense of inequalities 
(p9|)] resonator remains strictly discrete (i.e., having ex- 
actly zero line width) in the entire absence of dissipa- 
tion [sec Eq. ( |33bj )] despite the quite significant inter- 
mode scattering. This entirely correlates with a common 
knowledge regarding the spectra of conservative ballistic 
systems |^9|, So far as the wave equations for the 

resonator in question prove to be separable despite its 
boundaries' complex structure, one may (and should) re- 
gard such a resonator as the exactly integrable system. 
In Sec. [v], we provide the results of numerical calculations 
which enable one to make some definite conclusions of a 
general type about the chaotic properties of spectra of 
wave systems whose classical analogs are nonintegrablc. 



2. Strong scattering limit 

By strong intermode scattering we imply the limiting 
case pertinent to the large (as compared to unity) norm 
of the mode-mixing operator R in Eq. (|2^). Subject to 
strong mode mixing, each of the modes in this case be- 
comes quasiuniformly distributed over a large region in 
mode space. |5l| Fo r the resonator we study here this cor- 
responds (see Eq. ( |28b| )) to inequality 



—2 



(t/.s c ) 2 » 1 



(39) 



The asperities that obey this condition are referred to as 
"sharp" boundary inhomogeneities. 

In the limiting case (|9|), the operator potential ( p5| ) 
may be expanded in a power series of the inverse op- 
erator FT 1 . The expression located between projection 
operators in Eq. (E5[) is then transformed to 



(40) 



Here, the symbol Q^ v ^ denotes the diagonal (in mode rep- 
resentation) operator on M M whose matrix elements are 
trial Green's functions, Eq . (pfl|) , with all indices v ^ pi. 
Substituting expansion (Q) into potential ( p5| ) we ob- 
tain the following asymptotic expression for the diagonal 
mode propagator, 



G 



G 



(V) 



-(r~ 

2 V 



G 



(v) r 



1 



2 M 



' MM 



(41) 



The second term in square brackets in Eq. (|4l]) is a small 
(in the parameter ||R|| _1 <C 1) correction to the first one. 
Moreover, this corrective summand turns out to be less 
the closer the external frequency oj is to the resonance 
frequency of the "unperturbed" /^th resonator harmonics 
[the latter corresponding to vanishing the real part of 

Gfj, }. This implies that strong intermode scattering 
caused by the resonator boundary roughness does not 
have a great impact upon the spectral energies' analytical 
structure against the weak scattering situation. In both 
limiting cases the resonance frequencies are specified to a 
good parametric accuracy by the poles of the trial Greens 
function G^P , i.e., they can be found from the equality 
k « x M . 

At first glance, the width of the resonance lines be- 
haves in a somewhat abnormal fashion. In the absence 
of random inhomogeneities on the resonator side bound- 
ary it is determined by dissipation loss only, being equal, 
by the order of magnitude, to the rate l/rd; s (in units 
of "energy" A; 2 ). On sharpening the asperities within the 
range S < 1 the resonance line Q factors decrease as their 
widths are defined by the sum of the dissipation rate and 
the rate of scattering-induced mode decohcrence. 
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On passing to the region S > 1, as the notions of ray 
dynamics suggest, an increase of chaotic nature should 
be observed, which inevitably changes over to the ergodic 
stage. In the mode language this would imply the rear- 
rangement of the excited mode energy between a large 
number of resonator modes and, as a consequence, the 
quasihomogencous infilling of the resonator volume by 
the electromagnetic field, irrespective of the field struc- 
ture of the mode supposedly excited by the source. 

As a matter of fact, the result of the wave equation so- 
lution is different. Specifically, the resonator spectrum 
is back to its original, formally undisturbed form(|lg|) 
when the asperities are gradually sharpened. As to the 
disorder-induced composite mode de cay r ate [the analog 
of the weak-scattering formula Eq. ( 33b )], with infinite 
sharpening of the asperities it progressively vanishes. In 
the limit of S — > oo the width of resonance lines returns 
to the initial state, being governed by dissipation loss 
only. 

Such peculiar behavior of resonance line width cannot 
be interpreted in terms of classical, i. e. trajectory-based, 
chaos theory as the classic analog of the system we study 
here is definitely the K system p|. A decrease in the 
scattering-related portion of the resonance line width, as 
the random component of the resonator boundary in- 
creases, may solely be explained by the interference of 
scattered modes, whose contribution is lowered (accord- 
ing to the random-phase summation principle) with an 
increase in the number of interfering modes. 



V. NUMERICAL ANALYSIS OF THE 
ROUGH-RESONATOR SPECTRUM 

Now we present the results of numeric analysis of spec- 
tral properties of the rough resonator in the study, using 
the analytical formulas obtained in the previous section. 
The subject of the analysis is the statistical distribution 
of nearest level spacings (NLS) conventionally defined 
through unfolded spacings s^, 



s t » (E t+1 - Ei) g(Ei) , 



(42) 



where g(Ei) is the mean density of states evaluated at 
energy E t . 

According to the long-standing conjecture of Berry and 
Tabor ]42|] , the spectrum of an integrable system (nor- 
mally referred to as a regular spectrum) is described 
by the NLS distribution following Poisson's law, p(s) = 
cxp(— s). This distribution is normally associated with 
the absence of level correlation and, as a consequence, 
with their clusterization near s = 0. If there exists some 
chaotic component in the spectrum, the levels become 
correlated to a certain extent and subjected to a kind 
of repulsion, which is manifested through the gradual 
transformation of p(s) from Poissonian- to Wigner-type 
distribution, p(s) = (n/2)s exp(— irs 2 /4). The latter dis- 
tribution is typical for systems with extremely chaotic 
classical dynamics El. 



One can infer from our results, Eqs. ( J33| ) and (41), 
which apply to weak and strong intermode scattering, 
respectively, that in the idealized dissipation-free reso- 
nance system the line width should be zero, whatever 
the strength of the seed-mode scattering induced by sur- 
face inhomogencities. Yet, the mere presence of this type 
of scattering, even though it is weak, results in substan- 
tial redistribution of resonance frequencies in comparison 
with the case where boundary asperities are entirely ab- 
sent. The shift of each resonance line from its starting 
position (at S = 0) is determined by the sum of seed 
energy (Q) and the correction to this energy resulting 
from composite mode scattering. In the smooth rough- 
ness domain corresponding to 5 2 < 1 the correction to 
energy, Eq. dis]), has the form of Eq. (33a). In the op- 
posite limiting case, where asperities are strongly sharp 
(S 2 ^> 1), we regard this correction as negligibly small, 
using expression ( |4l"| ) with no corrective term for the /xth 

mode trial Green's function, . 

Figure ||, shows NLS histograms for inhomogeneity- 
frce resonator side-walls (3 = 0) and for some instances 
where parameter S is of finite but small value, as com- 
pared to unity. It is evident that if the dissipation loss 
is not taken into account (upper plots in Fig. |^) then 
the level spacings are distributed close to Poisson law, 
however rough the resonator may be. Considering the 
dissipation and, consequently, the resonance line finite 
width leads to partial "Wigncrization" of the NLS distri- 
bution, specifically, to the appearance of a significant dip 
of the function p(s) in the region of s ps 0. This is due to 
the fact that, when counting the number of resonances 
in a given frequency window, only resonance frequencies 
spaced at distances larger than the width of the corre- 
sponding lines should be taken into account. All other 
resonances, which are positioned very closely with a given 
one, in our calculations were considered as merging into 
a single resonance peak. 

On gradual sharpening the asperities, specifically, in 
the domain where S> 1, the resonance frequency distri- 
bution may well be described by formula (|||) only. The 
broadening related to r ough ness- induced mode scatter- 
ing (the analog of Eq. (33b) for S <C 1) monotonically 
decreases with growing parameter S, and in the limit 
S —> oo the width of all resonance lines is again deter- 
mined by solely the dissipation loss. Figure |3j shows the 
dependence of several composite mode energies on the 
asperity sharpness. With increasing parameter 5 all the 
energies monotonically increase. The intervals between 
neighboring levels tend to increase as well. 

This particular fact allows one to neglect (under the 
condition of strong scattering, S 1) not only the 
scattering-induced decoherence but also the initial dis- 
sipative broadening, which was chosen to be small from 
the outset so as to resolve the bulk of the resonance lines. 
As can be seen from the graphs in Fig. ||, the composite 
modes "float up" in energy as the asperities get sharp- 
ened. In each prescribed frequency interval the number 
of resonances decreases with growing S, going to zero as 
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this parameter tends to infinity. 

In Fig. |], NLS histograms for very sharp boundary 
asperities are given, when composite mode scattering is 
particularly strong. One can see that on sharpening the 
boundary irregularities the intcrfrcquency interval distri- 
bution, being of Poisson nature at small S and in the 
absence of dissipation, transforms gradually to the bell- 
shaped curve with a maximum placed at s = 1. With 
H — > oo this distribution, allowing for its normalization, 
approaches the S function, namely, S(s — 1). Such a form 
is characteristic for NLS distributions of a fortiori inte- 
grable systems, in particular, of a set of harmonic oscil- 
lators with commensurable frequencies. Note that this 
particular system can always be specified with a single 
governing parameter; i. e., it can be thought of to be, in 
a way, one-dimensional. 

Such a behavior of NLS distribution for a randomly 
rough resonator, which is at first glance unusual, can be 
explained in the following way. As an example, consider a 
resonator with no dissipation and follow the transforma- 
tion of p(s) with an increase in the gradient parameter S. 
At the early stage, in the region of small (as compared to 
unity) values of this parameter, the distribution retains 
Poissonian form (upper plots in Fig. ||), which may be 
accepted as a signature of small interlevel correlations. 
In this case, the set of energy levels of weakly rough two- 
dimensional wave billiard (we consider the case of a single 
z mode, k = 1) in any reasonably large frequency interval 
is governed by two parameters (two mode-index compo- 



nents), as one can see from Eq. (18) 



With growing the parameter S, the NLS distribution 
gradually transforms from the distribution of the Poisson 
type, with the maximum at s = 0, to the bell-shaped dis- 
tribution whose maximum moves up to the point s = 1 
(see Fig. |Z]). In doing so, in view of the above dis- 
cussed "float-up" effect, the system of resonance levels 
in any prespecificd frequency interval approaches ever 
more the one-parameter set, which is characteristic of 
regular physical systems. Indeed, the effective (gradicnt- 
renormalized) azimuth indices entering eigenenergy ( |l8| ) 
for oscillations with different Z's simultaneously go to 
zero when S unlimitedly increases, and resonances tend 
to become differentiated (quantized) by the radial index 
only. If azimuth indices of all oscillations were really 
strictly the same we would get the 5-shapcd distribution, 
p(s) = S(s — 1), instead of actual peaked distributions 
with finite widths. 

The finite width of the actually obtained peaked dis- 
tributions is accounted for by the peculiar fact that at 
arbitrary large but finite values of the parameter S one 
can always indicate frequency intervals (maybe, of quite 
large frequency values) where the bulk of resonance levels 
having the same azimuth index n and distinguished by 
radial index only is supplemented with a certain number 
of levels (the lesser is the number the larger is S) with dif- 
ferent azimuth indices. In the classical phase space such 
admixing of "foreign" levels results in the deformation of 
energy contours out of the phase plane, which normally 



is associated with the onset of chaos. Meanwhile, in the 
wave mechanics we deal with in the present work the res- 
onance levels, even though intermixed, in the absence of 
dissipation remain perfectly stable (i.e., of zero width). 
This is the main reason why we cannot speak about the 
presence of quantum (or wave) chaos in the rough-side 
resonator unless it is dissipative. 

The analogous situation was already discussed in the 
paper by Berry and Tabor (see Ref. 42 ), where the au- 
thors, using the definitely intcgrable system of harmonic 
oscillators as an example, have put forward and substan- 
tiated the concept of level clustering in the vicinity of 
either s = or s = 1. Obviously, the distribution of 
single-oscillator energy levels, which are entirely corre- 
lated, obeys the law p(s) — 6(s — 1). It was exactly the 
presence of other oscillators with incommensurable fre- 
quencies (the condition preventing the set of energy levels 
from being one-parametric) that has led the authors of 
Ref. to the NLS distributions similar in form to those 
depicted in Fig. || for two largest values of the gradient 
parameter, namely, for S 2 = 150 and S 2 = 300. 



VI. DISCUSSION 

In summary, we have suggested a new calculation 
method for finding the spectra of resonators with nonho- 
mogeneous boundaries. The method is based on direct 
solution of wave equations rather than on the system 
symmetry properties, and it is tried out for a cylindri- 
cal cavity resonator with a randomly rough side bound- 
ary (randomly rough two-dimensional wave billiard) . Al- 
though in the context of the classical-billiard-system the- 
ory this resonator is obviously nonintegrable, being con- 
sidered from a position of wave/quantum theory its gov- 
erning equation is shown to be separable, even if the sys- 
tem boundary is of quite complicated form. The condi- 
tion mandatory for the corresponding multidimensional 
wave equation to be separable is the violation of either 
V or T symmetry of the wave operator. For the inho- 
mogencities of random nature V symmetry is manifestly 
broken. This made it possible for us to study both dissi- 
pative and nondissipative rough cavity resonators within 
the unified approach. 

Using exact mode separation, which is carried out with 
our method at an arbitrary level of the disorder at res- 
onator boundaries, it is shown that the quality factors of 
resonance lines can become finite only if there is dissipa- 
tion and/or radiation loss in the system. In the absence 
of dissipation, despite the boundary complexity, the os- 
cillation spectrum of randomly rough resonator appears 
to be regular and possesses no signs of chaos. 

We also have found that the intensity of intermode 
scattering, and thus the degree of entanglement of res- 
onator modes, is predominantly determined by the aver- 
age sharpness of surface asperities rather than by their 
mean height. In the absence of dissipation and under 
conditions of weak oscillation scattering produced by the 
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resonator wall roughness (i.e., in the case of smooth 
boundary asperities) the interfrequency interval distri- 
bution in the resonator spectrum is close to Poissonian 
distribution. The resonance lines, being strictly dis- 
crete when there is no dissipation, with the presence of 
the latter become manifestly widened. The lines posi- 
tioned tightly on the frequency axis begin to merge with 
one another and cannot thus be considered as indepen- 
dent ones. The distribution of intervals between neigh- 
boring resonances broadened concurrently by dissipation 
and roughness-induced scattering acquires a significant 
Wigner component. The numerical estimations we have 
carried out on the basis of our formulas for nonsharp 
asperities, Eqs. ([33]), are in good qualitative agreement 
with our recent experimental results reported in Ref. . 

Within the suggested approach we were also able to 
consider, besides the case of smoothly rough resonators, 
the resonators whose boundary inhomogeneities are ar- 
bitrarily sharp. The sharpening of the asperities re- 
sults in the increase in energies of the effective (com- 
posite) oscillation modes, which gradually leave upward 
of any frequency interval chosen for spectral measure- 
ments. The spectrum of the quasioptical resonator, be- 
ing initially dense within any sufficiently large frequency 
range, in view of the mode "floating-up" effect that 
arises due to boundary inhomogeneities is substantially 
rarefied, thereby approaching monofrequency, or even 
zero- frequency regime if the asperities become sufficiently 
sharp. 

The results we have obtained in this study are per- 
tinent to billiard-type wave systems, which are ballistic 
from the viewpoint of classical dynamics. The conclu- 
sions we arrived at differ drastically from those obtained 
previously for resonators of similar cylindric form, though 
randomly inhomogeneous in the bulk infill rather than on 
the surface. The wave transport in disordered resonators 



with continuously varying bulk permittivity is always of 
diffusion nature P-fjlj , [p2[| irrespective of whether the 
dissipation is taken into account or not. For this reason, 
the state of chaos in their spectrum, which is determined, 
according to the Chirikov criterium, by the resonance line 
width jlEj, ^6|, is always finite. Hence, the conclusion re- 
sults that the spectra of wave and quantum resonance 
systems with continuous inhomogeneities in their bulk 
are always to some extent chaotic, whereas the analo- 
gous systems of truly billiard configuration reveal chaotic 
properties only provided that there exist some dissipation 
channels. 
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Appendix A: Evaluation of the mode-mixing 
operator norm 

When estimating the norm of the mode-mixing oper- 
ator R through Eq. (|27|), for the trial Green's function 
' we apply expression (20) with the intramode poten- 



cy 

tial V u put equ al to zero based upon arguments presented 
in Sec. jlYAl . 



As to the intcrmode scattering, to estimate the contri- 
bution to the norm (|27]) of the height potential V^ h ) one 

has to evaluate the correlator (ujyj^Uu)}^. Substituting 

the potential into Eq. ( p^ ) in the asymptotic form 
(|26|) and performing the averaging we obtain 



W(n v - n,j,)5„„>5 qvqfM 6i v i^ . 



(Al) 



Then, upon inserting expression Eq. (Al) into Eq. ( |27| ) and taking account of the normalization J2 n W(n) = 1 we 
arrive at the following estimate for the height-potential contribution to the intermixing operator R, 

(A2) 



<llRWf)~(jy «i. 



The correlator (^,/^*^i'/j L ' ) ) , which arises in estimating the norm of the gradient-potential operator R( sl ), after 



averaging can be recast as follows 

2 



(A3) 



It is easy to verify, subject to the particular form of eigenfunction (14c), that the product of integrals in Eq. (A3) is 
appropriately estimated by the following expression: 



R d f R d 1 

dr'ir'-J^l — lr';^)* \ dr{r; l v \ — |r; Z„) 



7i 



(l»vl) 



(I"mD 



dr 



i? 4 (I«mI) , (|n„|) . (|n„|) 

7; +7, 7, 7+7,, 



(A4) 



7j 



X 



7/.., 
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For estimation purposes we set n M ~ kR in Eq. (A4) and then, based upon the relationship -\/|n|(|n| + 2) < 7^ < 
■\/2|n|(|n| + 1) known from the Bessel function theory |ft7| , we use the inequality \f^ n ^ > kR. Substituting the 



asymptotic expression 



7; 



(l«l) 



2 2 



(Z » |n|) 



(A5) 



into Eq. (A4) inst ead of the exact roots 7; , one can make sure that replacing the product of integrals in Eq. ( |A3| ) 
with the product (A4) not only enables one to make a correct estimate of each of the terms in the sum (^) but also 
to carry out twofold summation over radial mode indices. The estimation formula for the operator R( sl ) norm is then 
reduced to 



1 /cr\2 



whereupon we arrive at the outcome 



<l|R-HV(0 



(A6) 



(A7) 



To estimate the norm of the operator R( s2 ) originating from the square-gradient potential (13b)one has to evaluate 
the following correlator: 



(n„) p(n„/) 



Ql * R 4 1 



B) B 



B 



n 2 



(I«mI) 



2 _ K, 



J 2tt J 2tt \i? 4 



(A8) 



The integrals over p and can be easily calculated with the use of the Gaussian model for the random process £,(p), 
and as a consequence we have 

/ / T^F "' 7^p" eXP -n^)ip + i{n v -n il )p'\W 2 {p- p') 

=2 (^) ^, n „ W(k)W{k + n M - n v ) . (A9) 



k—— oc 



In addition, if one takes the function W(ip) in 
the form of a Gaussian exponential, then M^(/c) = 
(ip c /\/27r) exp(— k 2 tp 2 /2), and the double integral over az- 
imuth angles in Eq. ( |A9| ) arc estimated by the simple 
formula 



/dpdp' 



exp 



(A10) 

The integrals in Eq. (A8), which contain the Bessel 
functions, can be estimated by substituting the asymp- 
totic expressions valid for large arguments for those func- 
tions. The sums over radial mode indices are then easily 



calculated, and one obtains the following estimate for the 
operator R( s2 ) square norm: 



(l|R (s2) H 2 )-(^)% c max^exp 



(?V - n v ) 2 p 2 



(All) 

In so far as the angle correlation parameter is small in 
magnitude, t p c <C 1, one can approximately replace the 
sum in Eq. (All) with the integral, thereby obtaining 



<I|R (S2) IIV(£) «!• (A12) 
From this estimation it follows that the quadratic in the 
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Cif ) potential in Eq. (|6|), in contrast to the linear poten- substitution into all of the formulas of the mean value of 
tial V^ 1 ), is adequately taken into account through the this potential. 
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The idea of such a smoothing goes back to Migdal's 
work on resonance levels of deformed nuclei, see mono- 
graph || 

The electric field of these oscillations is polarized parallel 
to the resonator butt-end plates, whereby these oscilla- 
tions are more sensitive to fluctuations of the resonator 
side boundary than are the TM-type oscillations. 
This assumption is justified to the measure of inequality 

©■ 

In classical chaos theory the characteristic size of this 
reg ion is specified as the dynamic localization length |39l 



The classic analogs of such resonators are nonconserva- 
tive dynamic systems, where energy is not conserved in 
time. 
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FIG. 2: NLS distribution for different values of slope parameter H 2 in the small-slope domain. Histograms are obtained from 
Eqs. (^). In each diagram, Wigner and Poisson distributions (solid and dashed lines, respectively) are shown for reference 
purposes. 
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FIG. 3: The dependence of composite mode energy (|l8|) on the degree of sharpness of side boundary asperities. 




FIG. 4: NLS distributions for large values of the slope parameter H 2 . Gradual level clustering in the vicinity of s = 1 is clearly 
seen in the case where mean sharpness of the asperities measures extremely large values. 



